Time Dependent Quark Masses and Big Bang Nucleosynthesis Revisited 
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We reinvestigate the constraints from primordial nucleosynthesis on a possible time-dependent 
quark mass. The limits on such quark-mass variations are particularly sensitive to the adopted 
observational abundance constraints. Hence, in the present study we have considered updated light- 
element abundances and uncertainties deduced from observations. We also consider new nuclear 
reaction rates and an independent analysis of the influence of such quark-mass variations on the 
resonance properties of the important ''He(d, p)*He reaction. We find that the updated abundance 
and resonance constraints imply a narrower range on the possible quark-mass variations in the early 
universe. We also find that, contrary to previous investigations, the optimum concordance region 
reduces to a (95% C.L.) value of 
averaged quark mass. 



0.005 < 5mq/mq < 0.007 consistent with no variation in the 



PACS numbers: 26.35.-l-c, 95.30.Cq, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

Big bang nucleosynthesis (BBN) remains as one of the 
key constraints on the physics of the early universe. BBN 
occurred during the epoch from 1 — 10* sec into the 
big bang, and as such, is the only direct probe of new 
physics during the important radiation dominated epoch. 
Moreover, once the nuclear reaction rates have been spec- 
ified, the light-element abundances produced during the 
big bang only depend [l| upon a single parameter, the 
baryon-to-photon ratio rj, which is now constrained to 
high precision from the effects of baryons on the second 
and third acoustic peaks in the power spectrum of cos- 
mic microwave background fluctuations as deduced from 
data obtained with the Wilkinson Microwave Anisotropy 
Probe (WMAP) ^. 

In this context there has been considerable interest 0- 
Q in recent years in the use of BBN to constrain any 
possible time variation of fundamental physical constants 
in the early universe. One can only detect variations 
of dimensionless quantities which are independent of the 



units employed. Therefore, of particular interest for the 
present work is the variation in the dimensionless ratio 



Xq = niq/AQCD, where 



md)/2 is the aver- 
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aged quark mass and ^oc d is the scale of quantum chro- 
modynamics (defined [5| as the energy associated with 
the Landau pole in the logarithm of the running strong 
coupling constant). As in Q we note that nuclear pa- 
rameters (e.g. nucleon mass, reaction cross sections, etc.) 
constrain t^qcD- Hence, it is convenient to assume that 
Aqcd is constant and calculate the time variation de- 
pendence in terms of 8mq/mq = 5Xq/Xq. Aqcd can 
then be added at the end to reconstruct a dimensionless 
parameter. 

A time dependence of fundamental constants in an ex- 
panding universe can be a generic result Q of the- 
ories which attempt to unify gravity and other inter- 
actions. Such grand unified theories imply correlations 
among variations in all of the fundamental constants 0] . 
However, it can be argued @ that BBN is much more sen- 
sitive to variations in Xq than other physical parameters 
such as variations in the fine structure constant. Hence, 
Xq may be the best parameter with which to search for 
evidence of time variation of fundamental constants in 
the early universe. It is of particular note that a recent 
study [8[ found that an increase in the average quark 
mass by an amount Smq/rriq = 0.016 ± 0.005 provides a 
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better agreement between observed light-element primor- 
dial abundances than those predicted by the standard big 
bang. A similar conclusion was reached in Because 
of the importance of such evidence for a changing quark 
mass, it is important to carefully reexamine all aspects 
of the physics which has been used to place constraints 
upon this parameter from BBN. That is the purpose of 
the present work. 

In addition to the interest in evidence for a variation 
of the physical constants in the early universe, it has also 
been suggested Q that such variations may provide 
insight into a fundamental problem in BBN. There is an 
apparent discrepancy between the observed primordial 
abundance of ^Li and that inferred from BBN when the 
limits on the baryon-to-photon ratio rj from the WMAP 
analysis are adopted. Indeed, standard BBN with the 
WMAP value for r/ implies an abundance of ^Li that is 
a factor of 2.4 - 4.3 (4-5 cr C.L.) times higher [l^l than 
the obscrvationally deduced primordial abundance. A 
great many possible solutions to this lithium problem 
have been proposed, (e.g. and references therein). 
Nevertheless, in view of the importance of the possibility 
that depleted ^Li indicated a variation of fundamental 
parameters, it is worth carefully reinvestigating this hy- 
pothesis. 

In this paper, therefore, we carry out an indepen- 
dent evaluation of the effects on BBN from a variation 
in the parameter SXq/Xq. In this effort we are moti- 
vated by new detailed analyses P, [l^l of the uncertain- 
ties in the observed light-element abundance constraints. 
We also make an independent evaluation of the resonant 
'^He((i,p)''He reaction rate based upon both the forward 
and reverse reaction dependence on dmq/mq. We find 
that, although the uncertainty in the results increases 
due to variations in the resonance parameters, the re- 
vised abundance constraints narrow the range of possi- 
ble variations in the quark mass from BBN. This latter 
constraint dominates and decreases the optimum concor- 
dance region to a value of —0.005 < 6mq/mq < 0.007. 
Although some variation in the quark mass is not ruled 
out, the results of the present study are consistent with 
no variation in the averaged quark mass. 



II. MODEL 

The calculations performed here are based upon the 
standard BBN network code of [l^ , with a number of 
reactions for light nuclei {A < 10) updated based upon 
the latest JINA REACLIB Database VI. [11], and/or 
the NACRE rates 0. The rate of the ^}ie{^ile^yBc 
reaction was additionally replaced with that of [131 • We 
also adopt the neutron lifetime of 878.5 ± 0.7 stat ± O-'isys 
from [i^l. The usually adopted world average value of 
885.7 ± 0.8 s is based upon an average deduced by the 
Particle Data Group [l^ in 2006 and repeated in 2010. In 
recent work [isj . however, it has been demonstrated via 
Monte Carlo simulations of the experiments that some 



of the data used in that average suffer from a —6 s sys- 
tematic error. Moreover, the lower neutron lifetime is 
based upon a number of improvements [l^] in the mea- 
surement. This shorter lifetime also better satisfies the 
unitarity test of the Cabibbo-Kobayashi-Maskawa ma- 
trix j20l | , and it improves the general agreement between 
observed primordial abundances and BBN pll |. 

As described below, we make use of an analytic scal- 
ing of the resonant '^Iie{d,p)'^}ie reaction rate from [2^ 
because it explicitly includes the relevant Q value and 
resonance energy dependences. We have modified the 
thermonuclear reaction rates based upon the procedure 
outlined in 0, Q with a refinement for resonant reactions 
as described below. In particular, we make extensive use 
of the relation between nuclear parameters and the stan- 
dard model of particle physics as delineated in ^ . 

A. Binding Energy Dependence 

Changing Xq affects nuclear reaction rates in several 
ways. The strongest influence on BBN abundances, how- 
ever, simply comes from the dependence of reaction rates 
on nuclear binding energies. The main effect of nuclear 
binding energies enters through the sensitive dependence 
of reaction rates on the reaction Q values in the exit chan- 
nel. This is because the Q value affects the energy of the 
outgoing products for reactions with a positive Q value. 
Whenever possible in the calculations described below, 
we modify the Q-value directly, Q ^ Q + 6Q in appli- 
cable formulas for the reaction rates. Nevertheless, for 
some rates, the Q-value dependence is not apparent and 
it is necessary to specify explicitly the effects of varying 
the Q value for the relevant reactions. 

The Q-value effects on the reaction rates are quali- 
tatively described in Q, which we now summarize and 
slightly generalize. For electric dipole radiative capture 
the cross section for photons of polarization m and energy 
can be written |23| as 

where IQiml is the electric dipole moment and v is the 
relative velocity in the entrance channel. At low energies 
the dominant Q-value dependence simply derives from 
the exit channel i?^ term Hence, one can write 

a{E) ^E^^^{Q + Ef , (2) 

where E is the reaction center of mass energy and E-^ is 
the energy of the emitted gamma ray. 

For low-energy reactions in which there are two nu- 
clear species in the final state, the reaction cross sec- 
tion depends upon the resonance penetrabilities in both 
the entrance channel (for resonant reactions) and the 
exit-channels [l^]. In the exit-channel the Q-value de- 
pendence arises from two sources. For one, the cross 
section becomes proportional to the final state velocity. 
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V (X {Q + Ey^"^ of the two emitted particles. For the 
other, one must also take into account the penetrabili- 
ties of the two exit-channel particles if they both carry 
charge. In this case the penetrability also defends upon 
the reaction Q value. Hence, one can write 



a{E) ^ (g + £;)i/2e-V^^o/(Q+s) 



(3) 



where Eq is the Gamow energy characterizing the exit- 
channel penetrability. 



Ec 



with 



(4) 



(5) 



while Z3 and Z4 are the exit-channel charges. 

At the temperatures associated with BBN {kT < 100 
keV), one generally has E <^ Q. Therefore, one can 
expand Eq. ([3]) in Q to find the dependence of the cross 
section on changes in 5Q^ 



cr = CTo 




(6) 



Although the Gamow term in Eq. ^ is usually small, it 
can be important for some reactions, e.g. ^Be(n,p)^Li or 
3He(n,p)3H 0. 

An additional sensitivity on the reaction Q values oc- 
curs for the reverse photodisintcgration rates which are 
related to the forward rates via the equation of detailed 
balance: 



oc (o-u)twd X e 



'Q/kT 



(7) 



For the special case of the p(n, 7)0 reaction, there is 
also an additional sensitivity to the position of a virtual 
level with energy = 0.07 MeV. This leads to (H] 



(8) 



As in [8|, we denote the sensitivity of nuclear binding 
energies to the averaged light-current quark mass by 



K 



5E/E 
5mq/mq 



(9) 



The K values for light nuclei have been summarized 
in 0, H, As in Q we use the values given by the 
AV18-I-UIX nuclear two-body interaction, with hadron 
mass variations calculated in terms of the rriq using the 
Dyson- Schwinger equation calculation of [26| . 



B. Effect on Nuclear Resonances 

We will treat the resonance reactions of BBN slightly 
different than that of Q . There are two reactions in BBN 
which are dominated by narrow resonances. These are 
the ^He(d, p)^He and ■^H(d, n)'*He reactions. The cross 
sections for these reactions can be factored into the usual 
separation of a penetrability and the astrophysical S fac- 
tor: 



a{E) 



(10) 



where S {E) contains information on the nuclear matrix 
element. For a resonant reaction, the 5*- factor can be 
written 



S{E)^ 



P{E) 



{E- Erf 



(11) 



where Er is the resonance energy and is the total 
resonance width. Any residual energy dependence be- 
yond that of the entrance-channel resonance (e.g. the 
direct capture contribution, and the penetrability in the 
exit channel) is usually contained in a polynomial P{E) 
that is fit to the measured cross section. We adopt the 
unperturbed cross-section parameters for these reactions 
in [1^ This gives resonance properties of Er = 

0.183 MeV and F^ = 0.256 MeV for the ^}le[d,pfYLe 
reaction, while Er = 0.0482 MeV and F^ = 0.0806 MeV 
for the ^H(d, n)^He reaction. The ground-state binding 
energies were taken from the experimental data in [Ij. 

Changes to the ■^He(d,p)^He reaction will affect the 
abundances of the primordial abundances of "^He and ^Be. 
Changes in the '^H((i, 7i)'*He reaction similarly affect the 
abundances of "^H and ^Li. There are, however, newer 
formulations for these resonant rates given in UJ and 
p^ . These revised rates incorporate newer data and an 
improved R-matrix fit to the data. However, these newer 
rates are given in an analytic form that does not explicitly 
manifest the total dependence on reaction Q value and 
resonance energy Er- To account for this, we adopt new 
rates but correct them using a scaling based upon the 
discussion above and the analytic reaction rate given in 

Explicitly, in [23 , the thermonuclear reaction rate for 
a resonant reaction can be written 



Na{<7v) 



[ 27 EgV^' 
\ AkT J 

Seff{Eo) 



■ cxp 



1 + [(£^0 - Er)/iTr/2W 



(12) 



where C — Na{8/ fJ-T^y^'^ is a normalization. For the re- 
actions with a two-body charged final state of interest 
here, the dependence on the exit-channel penetrability 
[Eq. ([3])] is contained in a polynomial expansion of S{E) 
in E. The peak of the Gamow window in the Gaussian 
approximation is Eq = E]j^ [kT /2Y/^ , so that the pene- 
trability in the exit channel is contained in a polynomial 
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For the "^116(6?, p)** He reaction, the resonance is an excited 
state of the compound nucleus ^Li*. For the '^H((i, n)^He 
reaction the compound resonance is ^He* . The observed 
unperturbed resonance excitation energies for Eq. (jl4p 
are i;?^. = -9.76 MeV and iJ^Hc* = -10.66 MeV. 

The resonance energies are related to the excitation 
energy in the compound nucleus and the net binding en- 
ergies of the reactants, i.e., 



J^{d,n) _ jt;^-^^^ — Ef — Ed . 



(15) 
(16) 



Unperturbed [Sniq/mq = 0) reaction rates {Na^ctv 
^H(d,n)''He reaction as given in Q (Cybu04), [i 



FIG. 1: 

for the 
(REACIO), and 
Tg = r/(io^ K). 



The net shift in the resonance position due to variations 
in niq is then the combination of the shift of three ener- 
gies, 

SEl^'P^ = 5E.U' - SEs^, - SEd (17) 
= {K'.U' E'^u- - Ks^^Es^^ - KdEd) 



Mi (Desc04) as a function of temperature ^ith the K values defined by Eq. dH) and summarized in 



expansion of Seff{Eo). Hence, combining Eqs. ([6]) and 
(fT2|) . a general scaling for any resonant reactions with a 
charged two-body final state can be written 



Here we point out that there is a consistency check on 
the Sniq sensitivity of the forward '^lle{d,p)^}ie reaction 
from the reverse '*He(p, d)'^Hc reaction. For this case we 
have 




[NA{av)] = [NA{av)]^x 

• l-f[(go-g°)/(r./2)F 
1 + [{Eo - S,.)/(r,/2)]2 



(18) 



(13) 



where [iVyi(crw)]o is the unperturbed reaction, and E'^ is 
the unperturbed resonance energy. 

The analytic forms for the unperturbed reaction rates 
given in [l5|, [ll] and [l^] appear to be quite different. 
Nevertheless, in spite of the improvements in experimen- 
tal data and theoretical fit adopted, the unperturbed 
rates are surprisingly similar. Figure [T] shows a com- 
parison of the three different unperturbed resonant rates 
over the range of temperatures of Tg = T/{10^ K) = 0.1 
to 2. Although there is some slight difference in these 
rates at high temperature (Tg > 0.5), in the temperature 
range of most relevance to BBN (Tg < 1), the three rates 
are essentially indistinguishable. Hence, in this study we 
have confirmed that there is no significant uncertainty in 
the constraints on Smq/niq introduced by the more re- 
cently deduced formulations [l^ and [IB] for the resonant 
reaction rates. 

As noted above, the dependence of the resonance en- 
ergy Er on variations of SrUq have been evaluated in 
Refs. g, Q and are adopted in the present study. To 
determine the sensitivity of the BBN resonant reactions 
to changes in Sniq, one begins by expanding the reso- 
nance energy in terms of the unperturbed and perturbed 
energies, i.e.. 



Er 



E2 + 5Er 



(14) 



SE^J'-'^^ = SEsw - SE4^^ - SEp 

= (A'sLj.E'sLi. — X4Hc£'4Ho — KpEp) 



As a test on the robustness of the constraint on Smq/niq 
we also include the variations in the resonance energy 
based upon the parameters in this reverse channel. For 
this case we have i^s^. = —3.131, /\5Hc* = —2.867, and 
the constraint: SEr^'"^^ — SEi^* — (S&jjo = for = n 
or p. This implies that -ftTs^. = Ki^ic ^ (^-^Ho/^^yi* ) = 
-1.08 X (-28.30)/£;5^.. 



Primordial Light-Element Abundance 
Constraints 



Before summarizing the results of our BBN calcula- 
tion and the constraints on Smq/niq, it is useful to first 
summarize the light-element abundance constraints. In 
this task we are much aided by a recent thorough review 
by locco et al. [l|, and also new constraints on the pri- 
mordial helium abundance (l2j . Our adopted abundance 
constraints are then as follows. 



1. Deuterium 

Deuterium is best measured in the spectra of narrow- 
line Lyman-a absorption systems in the foreground of 
high rcdshift Quasars. Unfortunately, only 14 such sys- 
tems have been found [H. Taken altogether they ex- 
hibit an unexpectedly large dispersion. This suggests 
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that there could be unaccounted systematic errors. This 
enhanced error can be approximately accounted for by 
constructing the weighted mean and standard deviation 
directly from the data points. Based upon this, we 
adopt the conservative range for the primordial deu- 
terium abundance of D/H = {2.87^^:11) x 10"^ from Q- 
This implies a 2cr (95% C.L.) concordance region of 

2.49 X 10"^ < D/H < 3.31 x 10"^ . (19) 

We note, however, that if one restricts the data to the 
six well resolved systems for which there arc multiple 
Lyman-a lines p7l | , one slightly lowers the 2a deuterium 
constraint to 

2.44 X 10"^ < D/H < 3.22 X 10"^ . (20) 

Adopting this lower constraint would strengthen the con- 
clusions of this paper. Nevertheless, we adopt the more 
conservative constraint of Eq. (|19l) based [l| upon all of 
the data. 



in Ref. [l[ based upon an analysis [s^ that included 
new observations and photoionization models for metal- 
poor extragalactic HII regions, along with new atomic 
physics computations of the recombination coefficients 
for Hel and the collisional excitation of the HI Balmer 
lines. However, the extraction of the final helium abun- 
dance has been fraught with uncertainties due to corre- 
lations among errors in the neutral hydrogen determi- 
nation and the inferred helium abundance. In [l^ it 
was demonstrated that updated emissivities and the neu- 
tral hydrogen corrections generally increase the inferred 
abundance, while the correlated uncertainties increase 
the uncertainty in the final extracted helium abundance. 
Therefore, we adopt the value and uncertainty from [l^ 
of = 0.2561 ± 0.0108, which is in general agreement 
with the predicted value from standard BBN when rj is 
fixed from the WMAP analysis 0. Hence, for the Yp 
constraint we adopt 

0.245 <Yp< 0.267 . (22) 



2. ^He 

The abundance of '^He is best measured in Galactic 
HII regions by the 8.665 GHz hyperfine transition of 
'^He+ [2^. A plateau with a relatively large dispersion 
with respect to metallicity has been found at a level of 
^Hc/H=(1.9±0.6) X 10"^. It is not yet understood, how- 
ever, whether ^He has increased or decreased through the 
course of stellar and galactic evolution [2l,[33|- Whatever 
the case, however, the lack of observational evidence for 
the predicted galactic abundance gradient (3l| supports 
the notion that the cosmic average "^He abundance has 
not diminished from that produced in BBN by more than 
a factor of 2 due to processing in stars. This is contrary 
to the case of deuterium for which the observations and 
theoretical predictions are consistent with a net decrease 
since the time of the BBN epoch. Moreover, there are re- 
sults [13 from 3D modeling of the region above the core 
convective zone for intermediate-mass giants which sug- 
gest that in net, ^He is neither produced nor destroyed 
during stellar burning. Fortunately, one can avoid the 
ambiguity in galactic '^Hc production by making use of 
the fact that the sum of (D + ■^Hc) /H is largely unaffected 
by stellar processing. This leads to a best estimate [l| of 
•^He/H = (0.7 ± 0.5) x 10^^ which implies a reasonable 
2(7 upper limit of 

^He/H < 1.7 X 10"^ . (21) 

and a lower limit consistent with zero. 



4. ''Li 

The primordial abundance of ''Li is best determined 
from old metal-poor halo stars at temperatures corre- 
sponding to the Spite plateau (see [l| and references 
therein). There is, however, an uncertainty in this de- 
termination due to the fact that the surface lithium in 
these stars may have experienced gradual depletion due 
to mixing with the higher-temperature stellar interiors 
over the stellar lifetime. On the other hand, there are 
limits on the amount of such depletion that could occur 
since most lithium destruction mechanisms would im- 
ply a larger dispersion in abundances determined from 
stars of different masses, rotation rates, magnetic fields, 
etc., than that currently observed. In view of this un- 
certainty a reasonable upper limit on the ^Li abundance 
has been taken [l| to be 6.15 x 10^ which is based upon 
allowing for a possible depletion of up to a factor of ~ 5 
down to the present observationally determined value of 
7Li/H= (1.2315;^^) X 10"^" (95 % confidence limit) [H. 
A lower limit can be taken from the 2cr observational un- 
certainty in the presently observed value. We here adopt 
this constraint [l| which implies a ^Li/H range of 

0.91 X 10"^" < ■^Li/H < 6.15 X 10"^" . (23) 



III. RESULTS 



3. ^He 

The primordial ^He abundance, Yp is best determined 
from HII regions in metal-poor irregular galaxies extrap- 
olated to zero metallicity. A primordial helium abun- 
dance of Yp = 0.247 ± 0.002stat ± 0.004syst was deduced 



In Fig. [5]we show primordial abundances as a function 
of variations in the quark mass Siriq/mq for a fixed value 
of 77 = 6.23 X 10-1° deduced from the WMAP 7 year data, 
i.e. n^K^ = 0m258ta°oooll for model ACDM-l-SZ-hlens. 
Our calculation of standard BBN {Smq/nig — 0) predicts 
the following primordial abundances: D/H=2.593 x 10"^, 
'''He/H=1.007 X IQ-^, Yp = 0.2466, 6Li/H=1.190 x 10"", 
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and '^Li/H=5.017 x 10-^°. We note that our calcu- 
lated standard BBN abundances differ slightly from those 
calculated previously in Values and correspond- 

ing error bars from [lo| for respective light nuclei at 
5mq/mq = are shown as points on Fig. [2] The un- 
certainties in the standard BBN abundances were esti- 
mated in [Toj from a Monte Carlo simulation incorpo- 
rating the uncertainties in nuclear reaction rates. Our 
standard BBN ^He abundance is smaller than that ob- 
tained in [lo| by an amount which exceeds the estimated 
BBN uncertainty. This difference mainly derives from 
our adoption of the shorter neutron lifetime from [l8| . 

The blue solid lines on Fig. [5] arc for the case of no 
shifts in the resonance energies as was also considered 
in 0. For these cases we have -fCsLj. = —1.54 and 

ifsHc- = -1-44 with 5Ei'^'^^ = 0. AU three resonant 
reaction rate evaluations of [l^ [l^, [22j were utilized. 
As noted above, however, the results for the different 
rate evaluations [H, [H, [22| are nearly indistinguishable 
from those obtained using the earlier rate deduced in [2^ . 
Hence, previous studies were justified in using the rates 
from [22j . The changes due to implementing the different 
rates are smaller than the drawn lines and much smaller 
than the uncertainties due to the errors in the BBN re- 
action rates as noted on Fig. [2] 

The dashed lines on Fig. [2] show results when the res- 
onance energies are shifted by the same amount as the 
ground-state binding energy shift so that the excitation 
energy in the compound nucleus remains fixed. In these 
cases we have JTsLi* = —3.35 and Ks^ic' — —3.19 and 
SEnA* = 6E&Ag^,, where ^Ag,s. denotes the ground state 
of the A=5 nuclei. 

The dot-dashed lines correspond to an average value of 
the resonance sensitivity based upon resonance parame- 
ters for the forward direction. For these curves we have 
= -2.29 and Ks^^. = -2.21. The sohd green 
curves are new to the present study. They derive from 
considering the reverse reaction rates for determining the 
variation of the resonance energy with quark mass. For 
this case we have i^sLj* = —3.131, /^shc* = —2.867, and 
the constraint: SEr^''''^ = SEsj^, — bEiYia = implies 
K^A- = ^*He X (Shc/^^^a*) = -1-08 X (-28.30)/£;5^.. 
We note that all lines for ^He are indistinguishable from 
each other, and that the dashed line for D is indistin- 
guishable from the green line. 

It is straightforward to understand the shape of calcu- 
lated abundances vs Smq/niq. As Smq/mq increases, 5Q 
for the ^H(n,7)^H reaction decreases. The neutron pro- 
cessing in BBN is then delayed. This leads to somewhat 
inefficient ^Hc production and a smaller "^Hc abundance. 
The neutron abundance is then higher, so that the ulti- 
mate D abundance is higher. 

Thus, the ''He abundance decreases and the D abun- 
dance increases as dniq/mq increases primarily as a result 
of a decrease in the ^H(rt,7)^H rate. The D abundance 
is, however, also slightly affected in the same direction 
by the '^He((i,p)*He rate. The "^He abundance increases 
as \Smq/mq\ increases. Changes in the resonance ener- 
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FIG. 2: Calculated light-element abundances as a function 
of variations in the quark mass Smq/niq for a fixed 77 = 
6.23 X 10"^° from the WMAP 7 year analysis Q]. The blue 
solid lines on this figure are for the case of no shifts in the res- 
onance energies as in §1, but for the three resonant reaction 
rate evaluations [isl, [Tg [23|. The dashed line corresponds to 
the resonances being shifted the same energy as the ground 
state. The dot-dashed line corresponds to an averaged value 
of the resonance sensitivity in the forward direction. The solid 
green curve is new to the present study. It derives from con- 
sidering the reverse reaction for the determining the variation 
of the resonance energy with Sniq/mq. The red boxes show 
the allowed parameter regions for the case of the reverse reac- 
tion determined using our adopted observational constraints. 
Theoretical uncertainties in standard BBN [l^ are shown by 
error bars at Sniq/mq = as a guide. 



gies in the reactions ^H(d, n)^He and ^He((i, p)^He lead to 
smaller rates as \5mqlm.q\ increases. This is because the 
energy levels of the resonant states, i.e., ^He* and ^Li*, 
arc farther removed from those of the energy of the ini- 
tial entrance channels 0. Thus, as \5mq/mq\ increases, 
the abundances of '^H and ^He increase. The abundances 
of ''Li and ''Be therefore also increase since the yields of 
''Li [via 4He(3H,7)] and ^Be [via *He(3He,7)] follow the 
production of ^ = 3 nuclei. 

From Fig. [2] it is also apparent that the revised con- 
straints on ''Li and "^He adopted here do not confirm a 
concordance best fit for SrUq/mq = 0.016 ± 0.005 as de- 
duced in Q. Rather, the optimum concordance level 
is for much smaller values of bmqjmq consistent with 
Smq/mq = 0. The constraint on the ^He abundance 
[Eq. (j22p ] limits the quark-mass variation to an upper 
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limit of 



5mjmq < 0.007 (^He), 



(24) 



The constraint from the D abundance [Eq. (fT9|) ] corre- 
sponds to 



0.005 < Smq/mq < 0.026 (D) 



(25) 



Note, also, that adopting the more tight constraint 
[Eq. (|20p ] from the six systems with well resolved multi- 
ple Lyman-a lines [23| slightly reduces the D/H limits to 
-0.007 < Sniq/mq < 0.023. 

There is no limit on Smg/mg from the '^He constraint 
[Eq. (PT|) ] if we allow for the possibility that there is no 
corresponding resonance energy shift in the '^H(c?, n)^He 
and '^He(d,p)^He reactions. On the other hand, if a res- 
onance energy shift with Smq/niq is allowed, then the 
■^He constraint could lead to an upper limit as small as 
Smg/mq < 0.006. 

The constraint [Eq. (P^j from ^Li leads only to a lower 
limit if we conservatively adopt only the averaged value 
of the forward resonance sensitivity (dot-dashed line on 
Fig. [11): 



0.005 < 6mq/mg (^Li) 



(26) 



Combining the above limits (except ^He), we deduce con- 
servative concordance upper and lower limits for Smq/rriq 
of 

- 0.005 < Smq/niq < 0.007 (conservative). (27) 

If more precise dependences of the rates upon the 
quark mass could be determined theoretically, a stronger 
limit on the quark-mass variation is possible. For ex- 
ample, if the case of a resonance shift from the reverse 
reaction (solid green line) is considered we obtain the 
limits of 

Smq/niq < 0.007 (''He), (28) 



0.003 < Smq/ruq < 0.018 (D), (29) 



Sruq/niq < 0.006 (^He), (30) 

- 0.004 < Smq/mq < 0.008 (^Li). (31) 

The combined concordance limit, however, is not much 
different from the conservative one [Eq. ^7} ]. This is 
because the quark-mass upper limit is primarily con- 
strained by the ^He abundance, while the lower limit is 
constrained from the D (or ^Li) abundance. Neither of 
these limits are strongly affected by the resonance en- 
ergy shifts. Nevertheless, including the resonance energy 
shifts strengthens the concordance constraints. 

We note that since this paper was submitted, it has 
come to our attention that Bedaque, Luu, & Platter 



[35| independently estimated the relation between quark 
masses and the binding energies of light nuclides. They 
reported a constraint of —0.01 < Smq/niq < 0.007 con- 
sistent with our conclusions. Their analysis, however, 
differs from ours in several ways. For one, their con- 
straint is based only upon the ^He abundance. More- 
over, they adopted a primordial helium constraint of 
0.240 ^ Yp < 0.258 that is more stringent than the recent 
value [Eq. (22)] adopted in the present work. Since no 
reference was given, it is difficult to assess whether this 
more stringent value is justified. Another difference is 
that the authors did not treat variations in nuclear reac- 
tion rates except for the ^H(n,7)^H reaction. Also, they 
adopted a neutron lifetime of 885 s. Although no refer- 
ence was given, this value is consistent with the old value 
[l9| , and different from the value 18 1 adopted here based 
upon more recent analysis. Although it is encouraging 
that they have reached a similar conclusion regarding a 
small value for Smq/mq, we believe that the results re- 
ported here are based upon a more recent and thorough 
analysis of abundance constraints and reaction rates. 

We note on Fig. 2 that the *He and D abundances 
exhibit a nearly linear variation with Smq/mq. Since 
these two abundances mainly determine the concordance 
region, it is useful for future reference to present here 
analytic formulas for the dependence of the ^He and D 
abundances with quark mass. For the primordial helium 
abundance we find 



Yp = 0.247-F2.0xlO"''(Tn - 878.5 s)-0.26 



5m„ 



ma 



(32) 



where the dependence of Yp on the neutron life t„ is also 
included. Note that this result is independent of which 
parametrization for the resonance shift is employed. For 
the deuterium abundance, in the case where the reso- 
nance shift in the '^H((i, n)^He and ^He((i, p)'*He reactions 
is obtained from the reverse reactions (solid green lines 
in Fig. dl, we obtain 



D/H = 2.59 + 39.2 



5m, 



X 10 



-5 



(33) 



while using the resonant shifts from the forward reactions 
(blue dot-dashed line) we obtain 



D/H= 2.59 + 34.6 ^ x 10 

m, 



-5 



(34) 



If we ignore the resonance shifts (solid blue line) we ob- 
tain 



D/H= 2.60 + 28.3- 



X 10" 



(35) 



Finally, it is worth remarking on the abundance of ^Li 
in the present studies. This nuclide is currently of inter- 
est as there has been some suggestion in observations [s^ 
of metal-poor halo stars for the presence of a primordial 
^Li abundance that is ss 10^ times the standard BBN 



8 



prediction. Indeed, ^Li production in standard BBN oc- 
curs via the ■*He((i, 7)^Li reaction and one might expect 
some consequence of a varying quark mass on its abun- 
dance. Nevertheless, the ^Li abundance remains smah in 
all of the parameter space considered here. The energy of 
the resonance is rather high (0.712 MeV) relative to the 
''He-l-d entrance channel. Moreover, its width is small, 
r = 0.024 ± 0.002 MeV [S^. Increasing the ^Li abun- 
dance by a significant amount would require shifting this 
resonance by ~ 0.7 MeV. This would, however, imply 
an unrealistically large Smq/rriq. For example, if the K 
value for the ^Li resonance is the same as that for the 
^Li ground state, the shift in the ''lle(d, 7)®Li resonance 
energy is given by 0, 

nig 

= 9.85 MeV^ . (36) 
rriq 

Hence, 5Er = —0.7 MeV would correspond to 5mq/mq = 
—0.07 MeV. As can be seen in Fig. 2, such a large shift in 
the quark mass would result in inconsistent abundances 
for ''He and D. Since this large deviation seems to be al- 
ready excluded, the effects of a varying quark mass on the 
resonant '*He((i, 7)^Li reaction can be safely neglected. 

IV. CONCLUSIONS 

We have reinvestigated effects of a hypothetical time 
varying quark mass on the light elements produced 



during big bang nucleosynthesis. The limits on such 
quark-mass variations are particularly sensitive to the 
adopted observational abundance constraints. Hence, 
in the present study we have considered updated light- 
element abundances and uncertainties deduced from ob- 
servations. We also consider updated reaction rates and 
an independent parametrization of the variation in the 
resonance energy in the ^He((i,p)'^He reaction from the 
reverse *He(p, d)'^He reaction. We find that the revised 
primordial abundances and constraints imply that there 
is no concordance best fit for 5mq/mq 0.016 ±0.005 as 
deduced in Q. Rather, the optimum conservative con- 
cordance region is for much smaller values of —0.005 ^ 
dniq/mq < 0.007, which is consistent with Smq/rriq = 0. 
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